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Abstract 

A recursion formula is derived which allows to evaluate invariant integrals over the orthog¬ 
onal group O(N), where the integrand is an arbitrary finite monomial in the matrix elements 
of the group. The value of such an integral is expressible as a finite sum of partial fractions 
in N. The recursion formula largely extends presently available integration formulas for the 
orthogonal group. 

PACS: 02.20.-a, 05.40.-a, 05.45.Mt 

1 Introduction 

Integrals over the classical compact groups 1,2 are of interest in various fields, such as harmonic 
analysis 3 or random matrix theory 4 . In these applications the integrand is often a polynomial in 
the matrix elements of the group itself, ( i.e. of the true matrix representation of the group). Thus 
we have to integrate an arbitrary monomial of these matrix elements. Closed formulas are available 
only for very special cases 3,5- ' and even a new method by Prosen et al. 8 using computer algebras 
is practically limitted to low degrees. Though note that, for arbitrary monomials there is strong 
evidence, that the results are exact at least up to the next leading order in V -1 with respect to the 
approximation of the group integral by independent Gaussian distributed matrix elements. 9 

In the present paper we shall address the case of the orthogonal group O(N). First we rederive 
the well known one-vector formula. 4,6 In this context, the terms “i?-vector formula” or “I?-vector 
integral” refer to the case where the monomial in question contains only powers of matrix elements 
from R rows or R columns respectively. Next we derive a recursion formula that relates an I?-vector 
integral to a linear combination of (R— l)-vector integrals. This is the central result of the present 
paper. Together with the one-vector formula, it allows to calculate any integral over a monomial of 
finite degree in a finite number of steps. This result is then used, to obtain a closed expression for 
general two-vector integrals that is much simpler than the one known before. 6 Besides, the older 
formula contains mistakes which (to the best of my knowledge) had never been corrected in the 
literature. 

The paper is organized as follows: In Sec. |2| we describe the current approach to the problem. In 
addition, we introduce some compact non-standard notations, which help to keep the mathematical 
expressions manageable. Then the one-vector result of Ullah 6 is rederived, as it is the base for 
the recursion formula developed later on. In passing we obtain an equally simple formula for the 
corresponding one-vector integral over the unitary group. In Sec. we derive the general recursion 
formula. In Sec. |] some applications are presented. As an immediate consequence, we obtain a 
closed expression for the two-vector integral, which is then compared to the corrected old result. 6 
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We also illustrate the use of our general formula for R > 2, calculating a particular three-vector 
integral. Sec. || contains the conclusions. 


2 General considerations 

To be specific, let us consider the orthogonal matrix w £ O(N) as a point in Euclidean TV 2 - 
dimensional space. Then we are interested in integrals of monomials in the coordinates of w. These 
are denoted by: 

N,R 

(M) = / dcr(w) H W^ . (1) 

J U=i 

Here er is the normalized Haar measure 10 of O(N), i.e. J d a(w) = 1, and M is a NxR, matrix 
of non-negative integers, with R < N. M is called the power matrix. In the recursion formula to 
be developed, R is used as the recursion parameter. Hence it is important, that i?, the number of 
columns of M, is as small as possible. 

The integral over the orthogonal group is invariant under any permutation of columns or rows of 
the integration variable w £ O(N). It is also invariant under taking the transpose. Therefore it is 
sufficient to consider such monomials which contain matrix elements from the first R < N columns 
of w only. According to Ullah 6 one may then write: 

< m )=-tttw > a/"(m )=[ n ( d °(^e) n^ i£ } n > ( 2 ) 

J ^ ’ ■' £=1 l *=1 J H<v 

where o is a NxR matrix of zeros. The integration region is the Cartesian product of R unit 
spheres with constant measures dfl(wj), and {wi,... , wr} are the corresponding unit vectors. The 
orthogonality of the unit vectors is implemented with the help of appropriately chosen 5-functions. 
{Wfi\w v ) denotes the scalar product between the two vectors: and w v . 


2.1 Compact notations for certain products of multinomials 


In the calculations to follow, we will frequently deal with certain products of binomial and multi¬ 
nomial coefficients. For better legibility we use two special notations: In what follows, x and y are 
TV-dimensional real vectors, and rh and n are TV-dinrensional vector-indices of non-negative integers. 
There are two typical cases in which products of multinomials appear: 

1) Consider the expression E = + Vi) ni ■ Its expansion gives: 


E = 



•''f' ki 


(3) 


where the sum runs over all k for which VT : ki < n t . In this case, the product of binomials is 
denoted by the following symbol: 


N 


n 




(4) 


2) In other occasions, we encounter expressions of the type: E' = where (r|u;') is 

the scalar product of the two (i?—1)-dimensional vectors r and w\. In this case the expansion 
reads: 

{ N 'j R -1 _ N 

nw*i>- • • * [ n r « 5 n w ^ ie ■ ( 5 ) 

t=l J c=l i= 1 


2 



Here we have to use the Nx(R-l) matrix K as an index. The elements of K are non¬ 
negative integers, k £ is the £th column vector of K , and k £ is the sum of its components: 

= Y^iLi The sum in Eq. (Q) runs over all I\ for which V* : -Kif = n,;. In this 

case, we use the following notation: 

N 

■ ■ ■ ,K itR - 1 ) = (n \K) . (6) 

2=1 


2.2 The one-vector integral 

In the one-vector case, R = 1, there are no orthogonality relations to respect. The power matrix 
M consists of one single column vector, here denoted by to. According to Eq. ®) we may write: 


( to ) 


M{fh) 

A M 


N(rh) 


N 

dn(«on^r, 

2=1 


(7) 


where o is a TV-dimensional vector of zeros. Following the original calculation of Ullah, 4,6 we 
integrate over the full vector space and implement the normalization of the column vector with 
the help of a (5-function. This introduces the integration constant Ci(TV): 


J\f(m) = Ci(TV) 1 



<5(||£|| 2 -1) . 


The next step is to remove the S- function. Setting Xi = it,/y^, we obtain: 


( 8 ) 


Af(m) r ( JV +™ )/2 - 1 = ci(TV ) -1 



<KN| 2 -r) , 


(9) 


where m = m *- The (5-function can now be removed by multiplying the l.h.s. and the r.h.s. 

with e~ r and integrating on r from 0 to cxd: 


V .oo 

A f(m) T (^±2) ci (TV) = [I / du 

-_i J — OO 


N 


=n r <^) 


( 10 ) 


Solving this equation for Af(m), the ratio N(m)/N(o) can be calculated, which leads to the desired 
result: 


N 

(^>=(tC/2 n (3)^/2 • ( n ) 

2=1 

Here it is convenient to use the Pochhammer symbol (z) n = T(z + n)/r( 2 :). 11 Note that Eq. (0) 
implies, that the integral (to) vanishes if at least one component of to is odd. 

2.3 The one-vector integral over the unitary group 

It is natural to consider also integrals over the unitary group U(N). This is in general much more 
complicated because usually the monomials to integrate contain powers of the matrix elements 
and their complex conjugated counterparts. However in the one-vector case, the integral over the 
unitary group can be mapped on a corresponding integral over the orthogonal group, which leads 
again to a simple result (it seems that such a formula has never been published elsewhere). If more 
vectors are involved, R > 1, the orthogonality conditions destroy this simple correspondence. 

To obtain the desired expression for one-vector integrals, it is convenient to consider monomials 
in the real and imaginary parts of the complex unit vector w. They can be identified with the 
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coordinates in a 27V-dimensional Euclidean space, where the Haar measure reduces to the constant 
measure O 2 on the unit hypersphere. Denoting the one-vector integral of an arbitrary monomial 


by (m:n) = Jl i=1 x T i V? i i where Wi = Xi + iyi, we may write: 


( m:n) — 


A4 (to, n) 
M(o, o) 


N 

M(m,n)= / d0 2 ('tZ;) 


x TvT 


Wi = Xi + iyi 


( 12 ) 


Note the different notations: (rn:n) stands for the one-vector integral over the unitary group, while 
(to, ft) is used in Sec. [| for the two-vector integral over the orthogonal group. 

Equation shows, that we may express (to : n) as a one-vector integral over the orthogonal 
group 0(2N): (to :n) = (p), where p is the 2iV-dimensional concatenation of to and n. Then we 
may apply Eq. ©■ This leads to: 

N 

(TO:n) = (iV)-^ +R)/2 [](i) roi/2 (i) rii/2 . (13) 

2=1 


Again the integral (m:n) vanishes, if at least one component of rh or n is odd. 


3 The recursion formula 

The desired recursion formula shall express an arbitrary integral (Af), where M is a power matrix 
with R columns, as a linear combination of simpler integrals ( M'), where M' has only R— 1 columns. 
Starting from Eq. (j|2j) one may attack this problem head on, and separate the integration on the 
last unit vector wr from the remaining integral: 


R -1 


N 


' R -1 


v(M) = / m du(«? € ) n 


w, 


Mu 

it, 


n««W» ,w R -i]rh R ) . (14) 


£=1 




Here rh R is the last column vector of the power matrix M, and 

,. ( R 'i r- 1 

,w R -i;fh R ) = dtt(w) n^((^)). ( 15 ) 

J U=i J 4=1 

As shown below, the value of this integral can be expressed as a linear combination of monomials 
in the integration variables {wi£ |£ < R — 1}. If this is inserted back into Eq. it obviously 
leads to the desired recursion formula. 

To evaluate the integral (|l5| ) , the integration over the unit sphere is replaced by an integration 
over the full space M. N , implementing the normalization condition with the help of a 5-function. This 
introduces again the normalization cons tant ci(N) [cf. Eq. (|[)]. The 5-function is then removed 
again, using the same trick as in Sec. 2.2. After that, the remaining 5-functions (responsible for the 
orthogonality relations) are replaced by their respective Fourier representations. In this way one 
obtains: 

{ N .. 1 R- 1 

P / d Xi xf iR \ 5 (||£|| 2 - 1 ) JJ 5 ({w^x}) 

1=1 ’’ ) 4=1 


= ci(N)- 1 r ( 


N—R+mR+l 




R -1 


r,MiR Xj 


n 6 


< i =1 * 

= ci(TV) -1 T (A^A±™a±i) _1 f J d Xi xf iR 

2=1 J 


4=1 


(16) 
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where w[ stands for the the row-vector (wn,... , wi^r- i) T - The integrals on Xi are easily evaluated, 
leading to: 


•/(-•-) r ( 


V-fi+mR+1 


iR— 1 


) cr(JV) = 


7 r 
N 

n 


R-l 


T e -T,^ =1 (r\0 


M iR 

^ ^ Kii even 

«i=0 


MiR ^ /•/-*! -*/\\MiR — K,i -p 

K, 




(17) 


Expanding the IV-fold product into a sum over the vector-index k, = (k,\, ... , kjv), we obtain for 
the l.h.s 


L/i.s. = ^ 


A™# 

V « 



(18) 


Note that the sum runs only over such k, for which all components are even, and that for the product 
of binomials the abbreviation from Eq. (||) is used. A bar over vector quantities such as rfiR and R 
denotes the sum of all their components. The quadratic matrix A , with elements A^ = (w^\w u ), 
has dimension R — 1. 

Now, the key observation is the following: The orthogonality conditions implemented in the 
form of ^-functions in Eq. & select from the total integration region only a sub-manifold. There 
it holds that (w^Wv) = 5^, so that the matrix A may be replaced by the unit matrix. Then it is 
possible to integrate the r-integral. The expansion of the product of scalar products leads to: 



where K is a matrix index with R — 1 columns, as introduced in Eq. (jb[) together with the abbre¬ 
viation for the product of multinomials. The fcj’s are the sums over the components of the column 
vectors of K. The remaining integrals are easily evaluated, and one obtains: 


J(...) 


r l -R 


E 


m R 

K, 


N 


(-i) 


ci(v)r( iV - il + 7Ttfl + 1 ) 

\ k ) n r (^o n 


( m H - re )/2 


n r m 


R -1 


N 


Kii 

w a 


K 


f=l 


i=1 


( 20 ) 


Note that, as a consequence of the r^-integrals, the sum in the second line runs over such K only, 
for which all k\,... Rr-\ are even. Inserting this expression into the initial Eq. ( 0 ) we obtain: 


A f{M) = 


A-R 


E 


( rriR 


c 1 (N)T( n - r+ ™* +1 ) YV ' 

^2(m R -K\K) (n 1 ^^) 


(- 1 ) ( 


N 


m R —K )/2 


n r m 


< i= 1 


R-l 


N 


X 

K 

{ R-l 

n i«v)) 

M< v > 


n / do (^e) n 


W, 


Miz+Kit 


( 21 ) 


£=i 


e=i 


i=i 
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The integral over the normalized vectors w 1 ,... ,iu r_i can be identified with J\f(M^ R ~ 1 ' 1 + K) 
which is a (R— l)-vector integral. In this way, we obtain a recursion formula for N(M). For the 
normalization constant, we find: 

*'■> %, (J v)r(I # n) {n r d)}{n r (i)}^- ) )' m 

Thus we obtain for the R -vector integral (M), defined in Eq. (gj): 

) (-1 )&«-*)/* jfj(I) 

x^rhn-KlK) ||l(I) Ss/2 | (M( r ~V + K). 

This is the desired recursion formula and the main result of the present paper. As mentioned 
before it is understood, that the first sum runs over such k only for which all components are 
even, while the second runs over such K only for which all are even. Furthermore 

rhu = Y2?=i Mm, R = and stands for the matrix consisting of the first R— 1 

columns of the matrix M. 

In principle Eq. ( |23| ) allows to evaluate integrals of arbitrary monomials of finite degree. The re¬ 
sult is always expressible as a rational function of the dimension N, because the repeated expansion 
of Eq. (^3|) leads to nested sums of partial fractions in N. In this context it is useful to note, that 
only the prefactor of the r.h.s. depends explicitly on N. The formula ( |23| ) can be used conveniently 
if either R or the degree of the monomial to be integrated are small. Otherwise Eq. (|2^) may lead to 
very lengthy expressions. However, such expressions should still be manageable with an appropriate 
computer algebra system. This would allow for further systematic studies of this class of integrals. 

In contrast to what one might expect, the integral (M) does not necessarily vanish if the power 
matrix M has odd elements. It rather holds the following: If any sum over a row or column of M is 
odd, then ( M) = 0. Though this is in fact well known, 12 it is still instructive to see that it follows 
almost immediately from the recursion relation (f23|). 

To this end, permute columns and rows, and take the transpose if necessary, to transform M in 
such a way that its last column contains the row or column whose sum of components is odd. Then 
apply Eq. (|X]): The sum over k is restricted to such k which have only even components. Hence R 
is even. As friR is odd, and ^5 = — k, a t least one of the sums k\,... , Rr-i must be odd 

as well. Such a term vanishes, because of what is said below Eq. (^o|). This implies that all terms 
of the sum over K vanish likewise, so that the complete integral gives zero. 

4 Applications 

In random matrix theory, many matrix ensembles are based on the concept of orthogonal invari¬ 
ance. Physically this corresponds to a situation where the Hamiltonian for a spin-less quantum 
particle possesses an anti-unitary symmetry, e.g. time reversal invariance. The Gaussian and circu¬ 
lar orthogonal ensembles 4,13 ’ 14 are well known examples based on this concept. Other examples are 
the Poisson orthogonal ensemble, 15 or the recently introduced matrix ensembles for semi-separable 
systems. 8 In those cases where the orthogonal invariance applies directly to the Hamiltonian, the 
statistical properties of the eigenvectors are uniquely determined by the orthogonal group and its 
invariant measure. Therefore any correlators between the eigenvectors can be expressed and calcu¬ 
lated in terms of i?-vector integrals. 


<«)=( t P)iE l 


/ rriR 
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In what follows, we first apply our integration formula ( p3| ) to the two-vector case. In this way 
we obtain a closed expression for arbitrary two-vector integrals. Then we compare this result with 
the corrected formula of Ullah. 6 For illustration, we finally compute a simple three-vector integral, 
which can be evaluated with an independent method also. As it should be, we find the same answer 
with both methods. 


4.1 The two-vector integral 

Consider the arbitrary two-vector integral (M) = (m,n), where the first column vector of M is 
denoted by m and the second by n. In this case, Eq. (^3|) leads directly to the following expression: 

= )«/2E(a) (- 1 ) ( ” _S)/2 |ri(5) Ki /2}(3)( S - B )/2 (m + n-K). (24) 


The sum runs over such k only, for which all components are even. A bar over a vector quantity 
denotes, as before, the sum of all its components. This formula has already been used in Ref. 16 to 
calculate various two-vector integrals. The numerical tests performed in parallel confirm its validity. 
For later purpose we use Eq. (^) to evaluate the following simple integral: 


(l 

1 

1 

\ 2 

1 / 

( 2 \ 
2 

0 

0 

/ N-l 

2 \ 

0 

V 




W 


-1 

(N - 1)N(N + 2) ' 


(25) 


Note that the same result can be obtained by an indirect method 12 also. 


In principle, an integration formula for general two-vector integrals has already been published 
some time ago. 6 After the correction of two misprints, it reads: 


n 2m< 2m _ -JV+1 ? -2iV-l-4-2£ (mi+ru) T(N — 1) T [N — 1 + + n i)\ 

i i F[E i m i + (IV-l)/2] + -l)/2] 

,2mjv,2ni,... ,2 niy 

X E (-1) E ^ (26) 

fei,... ,In— 0,... ,0 

II (\T) ( 2 D + h + l)/2] r[mj + th - (h + h- l)/2] 

T[N/2 + E.-fe + ii)/2] T[A^/2 + Ei( TO * + Ui) — A ^i)/2] 


where Vi : ki + h must be even. Here the original notation of Ref. 6 is used. The corrections concern 
the first line, where the nominator has been multiplied with T [N — 1 + + rij)], and the sum 

over fci,... , k N , h,... , In, where the vector-indices must start with zeros instead of ones. Finally 
the notation is quite unfortunate, as it seems to prohibit monomials with odd powers, though there 
is no reason for it. Indeed, Eq. (|g|) holds in those cases as well. This can be checked, for instance, 
by computing the integral (25) with the help of Eq. (^6|), setting m\ = m 2 = n\ = 712 = 1/2. Using 
the notation adopted in the present work, Eq. (Eq) reads: 


(m, n) = 


(-^ l)(m+n)/2 


(^l / 2 m 


4/2 


E 

k,l 


m 


( —1) ! (k + l) (rh— k + n — l) . 


(27) 


If we compare the integration formulas (|24|) and (^l]), they differ considerably. It seems rather 
difficult to prove their equivalence directly. Note moreover, that our result is much simpler, because 
there the sum runs over a single vector-index only. 
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4.2 A simple three-vector integral 

The three-vector integral considered here, is chosen because of its simplicity and because it may 
be evaluated using an indirect method, which allows to crosscheck the result. We will compute the 
integral (M) with 


M = 


(2 0 0 \ 
0 2 0 
0 0 2 
0 0 0 


\ : 


(28) 
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Henceforth we will skip those parts of the column vectors which are zero anyway. Using our recursion 
formula (|2^) the three-vector integral (M) can be reduced to a linear combination of two-vector 
integrals, for which we already have a closed expression, namely Eq. (]24|). Thus, the evaluation of 
(M) needs only a few steps: 



TV - 2 \(N- 1)N(N + 2) (TV - 1)TV(TV + 2)(TV + 4) 
TV 2 + 3TV — 2 

(TV - 2)(TV - 1)TV(TV + 2)(TV + 4) ' 


(29) 


The result 
nominator 


is expressed as a rational function, where care has been taken, that nominator and de- 
have no more common factors. 


Alternatively we may compute (M ), starting from the following identity: 


Z>« E 



= i 


(30) 


which holds for an arbitrary w € O(TV). Now we expand the products on the l.h.s. and integrate 
on both sides over the group. This gives: 

TV(TV — 1) (TV — 2) ^ ^0 2 0j ^ + 3TV(TV - 1) ^2 0 j ^ + TV ^2 j ^ = 1 . (31) 

It allows to express the three-vector integral (M) as a linear combination of a one-vector and a 
two-vector integral. According to the Eqs. & and (HJ), these integrals are given by: 



TV(TV + 2)(TV + 4) 



TV + 3 


0 2// l)N{N + 2)(N + 4) ' 


( 32 ) 
















Thus we finally obtain: 


I ° (N + 2)(7V + 4) — 3(iV + 3) + 1 

0 0 2j/ (N-2)(N-l)N(N + 2)(N + 4,) 

N 2 + 3N — 2 

~ (N - 2)(N - 1)N(N + 2 )(N + 4) ' 


(33) 


As expected, the result coincides with the one above, i.e. Eq. ( pg|). Here the indirect method worked 
so well because we first wrote down the identities (^) and and then chose our particular 

example (M). However, if the value of a certain integral is needed, one would have to guess useful 
identities which allow to express the integral by a linear combination of simpler ones, a procedure 
wich is certainly very difficult. In contrast to that, the recursion formula ( flS[ ) always provides a 
well defined finite procedure, for the computation of any integral. 


5 Conclusions 

To summarize, we have derived a recursion formula, which expresses an arbitrary I?-vector integral 
over the orthogonal group as a linear combination of (R— l)-vector integrals. It allows to successively 
evaluate the group integral of any finite monomial in the matrix elements of the group. The 
simplicity of the result depends primarily on R , the number of column or row vectors involved, and 
only secondarily on the degree of the monomial in question. The result is always given as a finite 
sum of partial fractions in N. 

As an immediate consequence of the general result, we obtained a closed integration formula 
for arbitrary two-vector integrals, which is quite different and much simpler than the corrected, 
previously known result. 

In principle a similar recursion formula can be obtained for integrals over the unitary group 
also. To that end one should consider monomials in the real and imaginary parts of the matrix 
elements. Though the derivation along the lines of the orthogonal case is rather straight forward, 
the resulting expressions are much more involved. It seem that the simple result for the case R = 1 
is only an exception. More work is clearly necessary to clarify the situation in this case. 
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